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AERONAUTIC  SYMBOLS 

1.  FUNDAMENTAL  AND  DERIVED  UNITS 


W 

S 

m 

I 


5 

6 
b 
c 

A 

V 
2 
L 

V 
D, 
Di 
D, 
O 


Symbol 

Metric 

English  , 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

Length _ 

1 

t 

F 

meter _ 

m 

s 

kg 

foot  (or  mile) _ 

second  (or  hour) - - 

ft  (or  mi) 
sec  (or  hr) 
lb 

Force - - 

weight  of  1  kilogram - 

weight  of  1  pound _ _ 

P 

y 

■ 

horsepower _ 

mph 

fps 

Speed _ 

fkilbmeters  per  hour - 

\meter3  per  second - 

kph 

mps 

miles  per  hour. _ 

feet  per  second _ 

2.  GENERAL  SYMBOLS 

Weig]it=m^  j  **  Kinematic  viscosity 

Standard  acceleration  of  bavity=9.80665  m/s*  ’  p  Density  (mass  per  unit  volunae)  / 

or  32.1740  ft/sec*  Standard  density  of  dry  air,  0.12497  kg-m  *-s*  at  15  C 

■pj7  and  760  mm;  or  0.002378  lb-ft“*  sec* 

Mass=^  Specific  weight  of  “sta.ndard”  air,  1.2255  kg/m®  or 

Moment  of  inertia=mF.  (Indicate  pis  of  0.07651  Ib/cu  ft 
radius  of  gyration  i  by  proper  subscript.) 

Coefficient  of  viscosity 

8.  AERODYNAMIC  SYMBOLS 


Area 

Area  of  wing 
Gap 
Span 
Chord 

.6* 

Aspect  ratio, 

True  air  speed 

1  tt, 

Dynamic  pressure,  2^1^ 
lift,  absolute  coefficient 


Drag,  absolute  coefficient 
Profile  di’ag,  absolute  coefficient 


D, 


Induced  drag,  absolute  coefficient 


D, 


Parasite  drag,  absolute  coefficient  Cdji=^ 
Cross-wind  force,  absolute  coefficient  Gc=^ 


G 


it 

Q 

Q 

B 


a 

e 

a„ 

at 

tto 


Angle  of  setting  of  wings  (relative  to  thrust  line) 
Angle  of  stabiliaer  setting  (relative  to  thrust 
line) 

Resultant  moment 
Resultant  angular  velocity 


Reynolds  number,  p 


VI 


where  Z  is  a  linear  dimen- 


sion  (e.g.,  for  an  airfoil  of  1 .0  ft  chord,  100  mph, 
standard  pressure  at  15°  C,  the  corresponding 
Reynolds  number  is  935,400;  or  for  an  airfoil 
of  1.0  m  chord,  100  mps,  the  corresponding 
Reynolds  nmnber  is  6,865,000) 

Angle  of  attack 

Anglo  of  downwash 

Angle  of  attack,  infinite  aspect  ratio 

Angle  of  attack,  induced 

Angle  of  attack,  absolute  (measured  from  zero- 
lift  position) 

Fhght-path  angle 
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SUMMARY 

A  surface  integral  is  developed  which  yields  solutions  of  the 
linearized  partial  differential  equation  for  supersonic  flow. 
These  solutions  satisfy  boundary  conditions  arising  in  wing 
theory.  Particular  applications  of  this  general  method  are 
made,  using  acceleration  potentials,  to  flat  surfaces  and  to 
uniformly  loaded  lifting  surfaces.  Rectangular  and  trapezoidal 
plan  forms  are  considered  along  with  triangular  forms  adapt¬ 
able  to  swept-forward  and  swept-back  wings.  The  case  of  the 
triangular  plan  form  in  sideslip  is  also  included.  Emphasis  is 
placed  on  the  systematic  application  of  the  method  to  the  lifting 
surfaces  considered  and  on  the  possibility  of  further  application. 

INTRODUCTION 

The  increased  emphasis  on  extending  theoretical  knowledge 
in  supersonic  wing  analysis  has  led  to  a  systematic  investi¬ 
gation  of  the  various  mathematical  methods  available  for 
treating  the  basic  differential  equations.  In  the  present 
report  advantage  has  been  taken  of  the  direct  analogy  whicli 
exists  between  the  linearized  partial  differential  equation 
for  supersonic  flow  in  three  dimensions  and  the  two-dimen¬ 
sional  wave  equation  of  mathematical  physics.  As  a  result 
of  this  correspondence,  solutions  which  have  been  given  for 
the  wave  equation  are  shown  to  be  applicable  to  the  type  of 
boundary  condition  encountered  in  wing  problems.  The 
first  section  of  the  report  is  devoted  to  the  development  of 
the  solution  for  the  potential  of  the  supersonic  flow  field. 
The  application  of  this  expression  to  a  number  of  examples 
in  supersonic  lifting-surface  theory  illustrates  the  usefulness 
of  such  a  method  of  attack.  In  the  first  of  these  examples 
the  loadings  over  the  given  plan  forms  are  assumed  to  be 
uniform.  The  results  obtained  for  such  cases  appear  at 
first  to  be  somewhat  academic  since  undesirable  twist  and 
camber  occur  over  portions  of  the  resultant  surfaces.  From 
the  uniformly  loaded  surfaces,  however,  it  is  possible  to 
develop  surfaces  having  arbitrary  load  distributions.  Im¬ 
posing  the  condition  that  the  final  lifting  surface  shall  be  a 
flat  plate  leads  to  the  solution  of  an  integral  equation  in 
every  case  considered.  The  results  obtained,  for  some  of 
the  plan  forms  considered,  have  been  developed  elsewhere 
but  not  always  with  the  unification  of  method  attained  here. 
New  configurations  are  also  included  among  the  examples 
given.  The  methods  shown  are  applicable  to  a  large  class  of 
unsolved  problems  of  immediate  interest. 


LIST  OF  IMPORTANT  SYMBOLS 


a 

local  velocity  of  sound 

A 

aspect  ratio 

b 

span  of  wing 

c 

chord  of  wing 

Cr, 

lift  coefficient 

C{0) 

load  distribution  function 

c. 

constant  value  of  discontinuity  in  ip  over  uni¬ 

formly  loaded  lifting  surface 

E(u,  k) 

incomplete  elliptic  integral  of  second  kind  with 

argument  u  and  modulus  k 

E,E' 

complete  elliptic  integrals  of  second  kind  with 

modulus  k  and  fl—E,  respectively 

F(u,  k) 

incomplete  elliptic  integral  of  first  kind  with 

argument  u  and  modulus  k 

IE 

functions  introduced  in  equations  (89)  and  (90) 

K,K' 

complete  elliptic  integrals  of  first  kind  with 

modulus  k  and  f  l  —  E,  respectively 

L 

lift  of  wing 

A/o 

free-stream  Mach  number 

7l,,7i.2,  % 

direction  cosines  of  normals  to  surface  S 

Vl 

static  pressure  on  lower  side  of  lifting  surface 

Pu 

static  pressure  on  upper  side  of  lifting  surface 

P:  (A,  Y,  Z) 

point  at  which  value  of  n  is  to  be  determined 

<L 

free-stream  dynamic  pressure  PoVo^^ 

S 

surface  enclosing  volume  V 

So 

area  of  wing 

U,  V,  w 

perturbation  velocities  in  direction  of  A',  Y, 

and  Z  axes,  respectively 

V 

volume 

Vo 

free-stream  velocity 

^,y,  2 

Cartesian  coordinates 

X,Y,Z 

transformed  coordinates  (See  equation  (3).) 

cn  {u,  k)  I 
dn  {u,  k)  > 
sn  {u,  k) ) 

Jacobi’s  elliptic  functions  of  argument  u  and 
modulus  k 

a 

angle  of  attack,  radians 

T 

Mach  forecone  from  point  P:  {X,  Y,  Z) 

A 

semivertex  angle  of  triangular  wing 
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2 

Ap 

S 

V 

e 

e{u) 

K 

\ 


A 


V\,  Vo,  Vs 

U(u,  y) 


Po 

<j,  9. 


T 

$ 

'P 

<Pu 


<PI 


CO 


pressure  differential  (pi—pu) 
angle  measured  from  X  axis 
conical  flow  coordinate  (See  equations  (27) 
and  (30).) 
i8  tan  5 

Jacobi’s  theta  function 

cylinder  of  infinitesimal  radius  enclosing  a.\'is 
of  forecone  r 

surface  at  which  stream  enters  induced  field  of 
wing 

angle  of  sideslip 

Mach  angle  of  the  free  stream  (^ii= arc  sin 

direction  cosines  of  conormal  v  to  surface  S 
incomplete  elliptic  integral  of  third  kind  with 
argument  u,  parameter  7,  and  modulus  k 
density  in  the  free  stream 

variables  representing  either  the  acceleration 
potential,  the  velocity  potential,  or  any  of 
the  three  perturbation  velocity  components 
surface  on  which  boundary  conditions  are  given 
velocity  potential 
acceleration  potential 

value  of  acceleration  potential  on  upper  side 
of  lifting  surface 

value  of  acceleration  potential  on  lower  side  of 
lifting  surface 


THEORY 

LINEARIZATION  OF  DIFFERENTIAL  EQUATION  FOR  COMPRESSIBLE  FLOW 


The  quasi-linear  (i.  e.,  linear  in  the  derivates  of  highest 
order)  differential  equation  for  the  velocity-potential  in 
the  case  of  compressible  fluid  flow  in  three  dimensions,  is 
expressible  in  the  form 


2^0 


a~ 


-2$, 


-2$  5^’"' — 0 


(1) 


where  a  represents  the  local  velocity  of  sound  in  the  medium 
and  Cartesian  coordinates  are  used.  Under  tlic  assump¬ 
tions  of  small  perturbation  theory  (references  1  and  2). 
this  equation  is  modified  so  that  it  is  linear  in  form  and 
consequently  more  amenable  to  mathematical  analysis. 
Denoting  by  the  variable  9  either  the  acceleration  potential, 
the  velocity  potential,  or  any  of  the  three  perturbation 
velocity  components,  the  linearized  expression  for  equation 
(1)  is 

(1 — (2) 


where  is  the  Mach  number  of  the  free  stream  and  thus 
equal  to  the  ratio  of  free-stream  velocity  and  the  correspond¬ 
ing  speed  of  sound. 


By  means  of  the  aflBne  transformation 
X=xx 

F=V±(l-Mo2)y  - 
Z=V±(l-Mo^)0. 


(3) 


equation  (2)  can  be  put  into  standard  forms.  Thus,  when 
Mo<  1  the  plus  signs  are  chosen  in  the  radicals  of  equation  (3) 
and  equation  (2)  becomes 

9xX'\-9YY~\~9};z  =  (i  (4) 

while  for  the  minus  signs  are  used  and,  as  a  conse¬ 

quence, 

9xx — 9yy — Slzz=0  (5) 

For  the  case  of  subsonic  flow  {Mo<C.l)  the  linearized  equation 
is  thereby  reduced  to  the  well-known  Laplace  equation  in 
three  dimensions.  Similarly,  in  supersonic  flow  (Mo>l) 
equation  (2)  is  again  reduced  to  classical  type  with  the  re¬ 
placement  of  the  space  coordinate  X  by  a  time  variable  T 
to  give  the  two-dimensional  wave  equation  of  mathematical 
physics.  The  linearization  of  the  general  differential  equa¬ 
tion  for  compressible  fluid  flow  therefoi’e  makes  available,  in 
both  subsonic  and  supersonic  studies,  the  results  of  the  ex¬ 
tensive  work  carried  out  in  previous  research  on  problems 
related  to  equations  (4)  and  (5). 

APPLICATION  OF  GREEN’S  THEOREM  TO  LINEARIZED  COMPRESSIBLE 
FLOW  EQUATION 

Methods  of  solution  for  partial  differential  equations  of 
the  type  considered  here  may  be  classified  into  two  principal 
categories:  methods  which  express  the  solutions  in  terms  of 
orthogonal  functions  and  methods  which  are  based  on  the 
use  of  Green’s  theorem.  Volterra’s  solution,  discussions  of 
which  may  be  found  in  references  3,  4,  and  5,  applies  the 
latter  approach  to  the  two-dimensional  wave  equation  and, 
as  a  consequence,  his  results  may  be  adapted  to  the  study  of 
supersonic  flow  and  specific  solutions  of  equation  (5). 

If  the  functional  notation 

i(0)  =  9xx  flyp  9zz 

is  used,  the  analytic  form  of  Green’s  theorem  for  equation  (5), 
relating  a  volume  integral  over  the  region  U  to  a  surface 
integral  over  the  surface  S  enclosing  V,  may  be  written  in 
the  form 


—  il,L{(T)\dV= 


((TDn9~9Dn<r)dS 


where  tr,  fi  ai'e  any  two  functions  which,  together  with  their 
first  and  second  derivatives,  are  finite  and  single  valued 
throughout  the  region  considered,  and 


D„9=ni 


where  ni,  n2,  ns  are  direction  cosines  of  inward  normals  to  the 
surface  S. 

The  expression  for  D„9  is,  of  course,  a  directional  deriva¬ 
tive.  The  corresponding  term  appearing  in  Green’s  theorem 
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for  Laplace’s  differential  equation  (incompressible  fluid  flow) 
is  precisely  the  directional  derivative  along  the  normal  to  the 
surface  S.  The  analogy  between  tbe  two  terms  prompts  the 
introduction  of  the  so-called  conormal  to  S  with  direction 
cosines  vi,  v^,  defined  as 

v\=—ni,v2—7i2,Vi—ni 

The  geometrical  connection  between  the  normal  and  the  co¬ 
normal  is  indicated  in  figure  1 ;  the  angles  between  the  lines 
and  the  Y  and  Z  axes  remain  respectively  equal,  while  the 
angles  between  the  lines  and  the  X  axis  are  supplementary. 


z 


FKii’RK  1.— Tho  gooinotric  relations  between  normal  and  conorrnal  to  surface  S. 


It  follows,  in  particular,  that  if  the  surface  S  is  the  plane 
the  two  lines  are  coincident;  if  (S'  is  a  cone  with  semivertex 
angle  equal  to  45°  and  a.xis  parallel  to  the  X  a.xis,  the  co¬ 
normal  at  any  point  lies  long  the  surface  S. 

It  is  now  possible  to  write 


_  dn  I  dn  bQ  bo 

LlnO-Ki  ^+^2 


(6) 


VOLTERKA’S  METHOD  FOR  TWO-DIMENSIONAI.  WAVE  EQUATION 

Consider  now  a  surface  t  which,  for  the  purposes  of  this 
report,  may  be  thought  of  as  being  coincident  with  the  XY 
plane  and  parallel  to  the  air  flow  which  is  in  the  direction 
of  the  positive  X  axis.  Two  such  surfaces  are  represented 
by  the  darkened  areas  in  figures  2  (a)  and  2  (b) .  It  is  desired 
to  determine  the  value  of  0  at  the  point  P:{X,  Y,  Z)  from  a 
knowledge  of  the  boundary  conditions  given  on  t.  The 
solution  to  such  a  problem  is  immediately  suggested  by 
equation  (8)  since  that  equation  requires  only  the  knowledge 

of  B  and  ^  along  a  surface  enclosing  a  given  volume,  together 
ou 

with  the  knowledge  of  some  particular  solution  o-  to  the  wave 
equation  valid  everywhere  within  the  enclosed  volume. 
Further,  it  is  physically  evident  that  contributions  to  the 
value  of  fi  at  P  can  come  only  from  points  within  the  forecone 
with  vertex  at  P  and  also  within  the  envelope  of  the  after¬ 
cones  with  vertices  at  the  foremost  disturbance  points  of  t. 
Referring  to  figure  2  (a) ,  this  would  mean  the  volume  bounded 
by  the  foreconc  T  and  the  wedge  X  springing  from  the  leading 
edge  of  t;  and  in  figure  2(b),  the  volume  bounded  by  the 
forecone  F  and  the  aftercone  X  with  vertex  at  the  apex  of 
the  surface  t.  Since  for  the  boundary-value  problems  in¬ 
volved  the  surface  t  remains  in  the  XY  plane,  equation  (8) 
must  be  applied  to  all  three  surfaces  X,  F,  and  t. 


and  the  surface-volume  relation  becomes 

J  J  J [<rP(n)-QP(<r)]cfF=  J  [[.  ds  (7) 

If  n  and  c  are  chosen  so  as  to  satisfy  equation  (5)  through¬ 
out  the  region  V,  then  equation  (7)  reduces  to  the  form 

The  form  of  equation  (8)  is  a  direct  analogue  to  results  ob¬ 
tainable  for  functions  satisfying  Laplace’s  equation.  (See, 
e.  g.,  reference  6,  p.  46.)  The  use  of  the  conormal  produces 
this  symbolic  equivalence. 


r 


Figure  2.— Concluded.  (6)  Triangular  plan  form. 
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Since  there  is  no  way  of  determining  U  and  ^  along  F  the 

attempted  solution  will  be  especially  difficult  unless  the 
particular  solution  a  and  its  derivative  with  respect  to  the 
conormal  vanish  everywhere  on  F.  But  this  is  in  fact  the 
essential  part  of  Volterra’s  method  of  solution.  Thus  the 
proper  choice  of  a  is 


£r=arc  cosh 


^(Y~Y,y+{Z-z,Y 


(This  relation,  incidentally,  is  the  indefinite  integral  of  the 
fundamental  solution  representing  a  supersonic  source  in 
three  dimensions  [(X-X,y-{Y-  Y,y-  {Z-Z,y]-''\}  The 
value  of  O'  is  equal  to  zero  on  the  forecone  F  since  the  equa¬ 
tion  of  this  cone  is 

ix-x,y- {Y-Y,y- {z-z,y=o 

and  further,  since  the  conormal  is  always  directed  along  the 
forecone,  ^  is  the  gradient  of  <r  along  F  and  is  also  zero. 

Equation  (8)  provides  an  equality  for  the  distribution  of 

0  and  ^  over  X  and  t,  provided  12  and  <7  satisfy  equation  (5) 

throughout  the  enclosed  volume  mentioned.  However,  al¬ 
though  <r  satisfies  equation  (5)  everywhere  in  the  enclosed 
volume  opposite  t  from  P  (under  the  XY  plane  in  fig.  2), 
along  the  line  (E— Zi)^=0  (above  the  XY  plane 
in  fig  2)  O'  is  infinite  and  does  not  satisfy  the  assumptions 
made  in  establishing  Green’s  theorem.  If  this  line  is  ex¬ 
cluded,  however,  by  means  of  a  cylinder  k  of  radius  e,  with 
axis  lying  along  the  line  (E~Ei)^-l-(Z— Zi)^=0,  then 
equation  (8)  may  be  applied  to  the  region  outside  k  and  yet 
within  the  space  bounded  by  X,  t,  and  F.  In  fact  equation 
(8)  can  then  be  written 


where  tj  is  the  portion  of  r  bounding  the  region  of  integration. 
If  R=  V (.Y—  Fi)^-1-  (Z—Z,y  and  cylindrical  coordinates  «,  4', 
and  (X~Xi)  are  used,  an  element  of  area  on  the  cylinder  k  is 
dS= — ed4'd(X — Xi),  while 


d(r_  btr  (F — X\) 


so  that 


hm  f  dS=\xm  f 

.->oJ  JA  bvj 

limj"  J^e  ^  arc  cosh  d4/d{X—Xi)  = 

J  (in 

=  -2x  r^l2(f,F,Z)df  (1 

JX). 


If  this  result  is  applied  to  equation  (9),  one  gets 

is  (11) 

and,  after  differentiating  equation  (11)  with  respect  to  X, 

PROCEDURE  FOR  LIFTING  SURFACES  AND  SYMMETRIC  WINGS 

When  the  region  considered  is  that  bounded  by  the  surface 
T,  F,  and  X',  the  portion  of  X  on  the  opposite  side  of  t  from 
the  point  P,  then  a  is  finite  throughout  the  region  and,  as  a 
direct  consequence  of  equation  (9), 

where  Q'  is  the  value  of  the  potential  function  on  the  side  of 
r  opposite  P  and  v'  is  in  the  opposite  direction  to  y  on  r. 
Adding  equations  (12)  and  (13), 


Q(X,Y,Z)-= 


1  d 

2ir  dA 

A  A 

2r  dX 


512  .  do' 


adS-{- 


The  integrations  over  r  are  now  in  a  form  which  may  be 
interpreted  directly  in  terms  of  known  conditions  over  bodies 
with  given  load  or  symmetrical  section.  The  integration 
over  X  and  X'  can  be  disposed  of  by  discussing  the  two  cases 
shown  in  figure  2.  When  12  is  identified  with  the  velocity 
potential,  its  value  can  be  shown  to  bo  zero  on  X  and  X' 
regardless  of  whether  the  leading  edge  is  swept  ahead  of  or 
behind  the  Mach  cone.  Wlien  12  represents  acceleration 
potential  or  any  of  the  perturbation  velocity  components,  a 
discontinuity  exists  in  the  value  of  12  for  leading  edges  swept 
ahead  of  the  Mach  cone  as  in  figure  2  (a) .  Analysis  of  this 
case,  however,  reveals  that  for  all  wing  problems  the  inte¬ 
gration  over  X  just  cancels  the  integration  over  X'.  When 
the  leading  edge  is  swept  behind  the  Mach  cone  as  in  figure 
2  (b)  the  value  of  0  is  again  zero.  Thus  in  any  case  there 
results  the  fundamental  equation: 


Q(X,Y,Z)- 


■  zA  A 

277  bA 


(12-12')  ^  dS 


The  counterpart  of  equation  (14)  for  incompressible  fluid  flow 
is  well  known.  (See,  e.  g.,  p.  60  reference  6.) 

Under  the  particular  conditions  for  which 
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over  the  surface  t  equation  (14)  becomes 

(16) 

The  restrictions  imposed  in  equation  (15)  can  be  given 
physical  significance  after  the  functions  Q,  O',  and  the  surface 
T  have  been  given  specific  meanings.  Consider  first  the 
case  where  r  is  a  lifting  surface.  Obviously  the  normal 
induced  velocity  w  is  a  continuous  function  across  t.  If 
0  and  O'  are  velocity  potentials  associated  with  the  lifting 
surface, 

ip(;?,  r,  z)  (X,  Y,  2) =5^= 

and  equation  (15)  is  satisfied.  If  0  denotes  acceleration  po¬ 
tential  or  perturbation  velocity  u,  it  is  necessary  to  show  that 
on  the  lifting  surface 

du_  du' 

This  relation  holds,  however,  for  since  w(X,  Y,  Z)  =w'{X,  Y,  Z) 
along  T,  it  follows  that 

and  from  the  condition  of  irrotationality  it  is  possible  to  ex¬ 
press  the  gradient  of  w  in  the  X  direction  as  the  gradient  of 
u  normal  to  the  surface,  that  is,  in  the  directions  of  v  and  v' . 
Equation  (16)  is  thus  applicable  directly  to  lifting-surface 
theory  in  conjunction  with  either  velocity  or  acceleration 
potentials.  Application  can  also  be  made  to  the  determina¬ 
tion  of  pressure  distribution  over  the  surface  of  a  symmetric 
airfoil  at  zero  angle  of  attack.  In  this  so-called  nonlifting 
case  the  function  is  set  equal  to  the  induced  velocity  w, 
r  is  the  plane  of  symmetry  of  the  airfoil,  and  equation  (16) 
can  be  used  to  establish  the  boundary  conditions,  provided 
equation  (15)  is  satisfied.  For  this  to  be  so  bwibv  must 
equal  —bw'jbv'.  But  conditions  of  symmetry  give  w(Z)  = 
—  w'( — Z)  from  which  the  equality  is  seen  to  hold. 

RETRANSFORMATION  OF  COORDINATES 

Since 

direct  substitution  into  equation  (16)  yields 
n(X,F,Z)  = 


a(x,y,z)  = 

^  p  r  _ (12—0')  jx—Xj)  {z—Zi)dX]dyi _ 

2jr2»zJ  Jr,  [(j,_y,)24-(2_2,)2]7(a;— a;j)2— /32[(i/— i/,)2-l-(3— 2i)2] 

(17) 

where 

/32=Mo'‘-1 

APPLICATIONS 

GENERAL  REMARKS 

Applications  in  lifting-surface  theory  may  proceed  along 
two  possible  lines  depending  upon  the  boundary  conditions 
specified.  In  what  is  usually  referred  to  as  the  direct 
problem,  or  problem  of  the  first  kind,  the  loading  is  given 
over  the  wing  and  the  potential  function  of  the  flow  field 
field  is  calculated.  From  the  potential  function  the  shape 
of  the  aerodynamic  surface  supporting  this  load  can  be 
found  relatively  easily.  The  inverse  problem,  or  problem 
of  the  second  kind,  concerns  itself  with  the  determination  of 
the  loading  over  a  wing  surface  from  a  knowledge  of  the 
surface  shape.  In  the  following  sections  both  of  these  cases 
will  be  considered.  The  direct  problem  will  bo  discussed 
for  various  plan  forms,  the  analysis  proceeding  directly  from 
the  expression  for  the  potential  function  given  in  equation 
(17).  The  detailed  discussion  of  the  direct  problem  is  justi¬ 
fied  by  its  application  to  the  inverse  problem  where  the  load¬ 
ing  over  flat  plates  with  rectangular,  trapezoidal,  and  tri¬ 
angular  plan  forms  is  determined.  The  mathematics  of  the 
inverse  problem  is  less  straightforward  since  the  analysis 
involves  the  introduction  of  elemental  lifting  surfaces  with 
constant  loading  and  the  solution  of  an  integral  equation 
for  each  plan  form. 

UNIFORMLY  LOADED  LIFTING  SURFACES  IN  SUPERSONIC  FLOW 

Infinite  span  wing, — In  order  to  determine  the  induced 
velocities  on  the  surface  of  an  infinite  span,  uniformly  loaded, 
supersonic  lifting  surface  by  means  of  the  methods  derived 
in  the  preceding  section,  it  is  convenient  to  set  O  equal  to 
the  acceleration  potential  tp  (reference  2) .  The  lifting  surface 
is,  in  this  case,  a  surface  of  discontinuity  for  the  function  tp 
and  corresponds  to  the  surface  ti  in  equation  (17).  The 
discontinuity  in  the  value  of  <p  between  the  upper  and  lower 
surface  is  equal  to 

(Fi— Z’u) 

Po 

where 

Po  density  in  the  free  stream 
Pt  static  pressure  on  lower  surface 
Pu  static  pressure  on  upper  surface 


1  f  f  (n-n')  (X-X,)  {Z-ZMXdY, 

2ir  bXj  j.,  [(F- Yif+  iZ-ZY] ^{X-X,f-(Y- Y,Y- (Z-Zyf 

This  solution  applies  to  equation  (5)  and,  in  order  to  relate 
problems  to  the  linearized  equation  (2),  it  is  necessary  to 
use  the  transformation  of  equations  (3).  If  the  point  X^, 
Fi,  Zi  transforms  to  the  point  Xi,  yi,  Zi,  it  follows  that 


It  follows  that  for  the  uniformly  loaded  wing  in  the  plane 
2i=0  the  discontinuity  in  the  acceleration  potential  is  a 
constant,  say  Go-  From  equation  (17) 


¥’(a:,2/,2)  = 


^  ^  r  r _ {x—xi)zdxidyi _ 

2x  bxj  j  [(^—y^y^z^]^/{x—Xly—^^[(y—yly+z^\ 


(18) 
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Mo 


Figure  3.— Regions  of  integration  for  infinite  span  unswept  wing. 

A  sketch  of  the  airfoil  plan  form  is  given  in  figure  3  and 
two  possible  regions  of  integration  are  indicated.  In  all 
cases  the  integration  with  respect  to  y  is  performed  between 
the  limits  at  which  the  radical 


■^1  {x—XiY— +  z^] 


vanishes  while  the  integration  with  respect  to  x  depends  upon 
the  manner  in  which  the  forcconc  of  the  point  P  intersects 
the  discontinuity  surface.  Denoting  the  chord  length  of  the 
airfoil  by  c,  the  following  relations  are  obtained: 

^=0  when  a;T/32<0 

v5=  (7o  when  (19) 

^  =  0  when  c<Yx^P2 

(When  double  signs  are  used,  the  upper  sign  refers  always 
to  the  case  where  2>0  and  the  lower  sign  corresponds  to 
s<0.) 

The  value  of  the  acceleration  potential  is  thus  seen  to  be 
zero  at  all  points  in  space  except  for  those  points  lying  within 
the  region  between  the  wedges  extending  back  from  the 
leading  and  trailing  edges  of  the  airfoil. 

It  is  now  possible  to  determine  the  induced  velocities 
associated  with  the  acceleration  potential  just  obtained. 
Since,  in  linear  perturbation  theory  (reference  2), 


dy~^°dx’  dz~'^^  dx 


(20) 


where  u,  v,  w  are  respectively  the  x,  y,  z  components  of  the 
perturbation  velocities,  it  follows  that 

d  f  ''  1 

Vo 


The  induced  velocities  for  the  infinite  span  airfoil  result 
immediately  from  equations  (19)  and  (21).  If  the  upper 
sign  of  a  double  sign  is  again  referred  to  the  2)>0  case,  the 
results  may  be  written  in  the  form 


v=0 


Mor  0:^X^Pz<^C 


(22) 


Since  the  vertical  induced  velocities  are  constant,  it  follows 
that  the  supersonic  airfoil  of  infinite  aspect  ratio  and  uni¬ 
form  load  distribution  is  a  flat  plate.  The  relations  between 
this  loading  and  angle  of  attack  will  be  considered  later. 

Lifting  surface  with  rectangular  plan  form. — The  complete 
discussion  of  the  supersonic  lifting  surface  with  uniform 
loading  and  rectangular  plan  form  is  lengthened  considerably 
by  the  fact  that  in  calculating  the  acceleration  potential  at 
the  point  P  with  coordinates  x,  y,  z  it  is  necessary  to  distin¬ 
guish  between  several  regions  in  space  in  which  the  point 
may  be  located.  These  regions  arise  from  consideration  of 
the  manner  in  which  the  forecone  of  the  point  P  cuts  the 
surface  of  discontinuity.  The  value  of  <p  can  be  found  with 
approximately  equal  facility  in  each  of  these  regions  but, 
since  this  paper  is  concerned  primarily  with  effects  on  the 
surface  of  the  airfoil,  the  solutions  for  pertinent  regions  only 
will  be  given  here. 

Figure  4  shows  the  rectangular  plan  form  LL'  T'  T  together 
with  the  coordinate  system  to  be  used.  The  dimensions  of 
the  wing  are  chosen  so  that  the  Mach  cones  extending  back 
from  the  leading  edge  will  not  intersect  within  the  boundaries 
of  the  wing.  This  restriction,  which  is  not  necessary  but 
merely  simplifies  the  analysis,  implies  that  if  b  is  the  span 
of  the  wing  and  c  the  chord  length,  then 

(23) 


where  /x=arc  sin  is  the  so-called  Alach  angle  of  the 

stream  and  equal  to  the  semiverlex  angles  of  the  Mach  cones. 
The  loading  over  the  rectangular  plan  form  is  to  be  uniform 

so  the  expression  —  tpi  is  s(‘l  equal  to  t  for  — 2  — 2  ^ 

and  0<Cx,Vc.  The  acceleration  potential,  expressed  as  a 
function  of  x,  y,  2,  is  thus  obtainable  from  equation  (17)  and 
the  limits  of  integration  must  be  determined  from  the  position 
of  P.  From  reasons  of  symmetry,  only  the  portion  of  space 
for  which  yVO  need  be  considered.  Once  the  acceleration 
potential  has  been  calculated,  equations  (21)  may  be  used 
to  calculate  induced  velocities.  The  results  of  such  calcula¬ 
tions  are  given  and  the  same  convention  for  double  signs  is 
used. 

Region  I, ;  Behind  the  leading-edge  wedge,  ahead  of  the 
trailing-edge  wedge,  and  bounded  laterally  by  the  y~0 
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Fiqukb  4.— Lifting  surface  with  Mach  cones  and  coordinate  system  for  rectangular  plan  form 

plane  and  the  Mach  cone  from  the  leading-edge  tip.  The 
results  in  this  region  correspond  to  results  obtained  for  the 
infinite  span  airfoil.  Thus 

<P=  ±  2  ^0 

u=  ±  6o/2T^o 
e;=0 


(24) 

Region  R:  Within  the  Mach  cone  from  the  leading-edge 
tip,  outside  the  Mach  cone  from  the  trailing-edge  tip,  and 
forward  of  the  trailing-edge  wedge.  Denoting  the  integrand 
in  equation  (17)  by  the  symbol  /,  the  expression  for  tp,  when 


(^40 


v= 


2,rF„ 


—z 


{y-h 


+  z^ 


Co  . 
'*®“27rFc 


—  arc  tan 


+  (26) 

{y~\  +2" 


As  a  partial  check  of  the  expression  for  (p  in  equations  (26) , 
it  can  be  seen  that  in  the  limit  as  z  approaches  zero  the  value 
of  ip  agrees  with  the  result  given  in  equation  (24)  on  the 
wing  while  the  value  is  zero  off  the  wing. 

The  values  of  vertical  induced  velocity  in  the  plane  2=0 
are  of  particular  interest  since  from  a  Itnowledge  of  the 
distribution  of  w  the  surface  shape  and  local  angle  of  attack 
corresponding  to  the  imposed  load  distribution  can  be 
determined.  The  expressions  for  w  for  uniform  loading 
wiU  be  particularly  useful  later  when  the  load  distribution  is 
modified  in  order  to  obtain  airfoils  with  specified  induced 
velocities.  Introducing  the  notation 


V 


(27) 


y<2  is 


'”(’=-»'")  =  +2Va(r'*-Jn  Jr?-'*') 

where 

Yi=y—^  ■yIix—Xiy—0^z‘^,  Xi=x—p^(y—^  +z 

^2=2/+w  V(a:— a;i)^— /3V,  X2=xYfiz 


(25) 


Application  of  equations  (21),  after  integrating  either  equa¬ 
tion  (25)  or  its  companion  expression  when  y'^2  yi®ids  the 
results: 


-arc  tan 


x{ 

y-\h) 

•* 

(y-\  b 

the  following  results  are  obtained  for  the  area  covered  by 
the  tip  cone: 


For  0<ij<(l, 


-GoP  n-y/l-v? 
2irFoJi 


dri 


1 


and  for  — 1<7;<0,  (x<c) 


^2=0 


-CoP 

2xFo 


(28) 


After  integration  of  these  two  expressions,  the  explicit  value 
of  vertical  induced  velocity  throughout  the  entire  region  is 
found  to  be 


Equations  (28)  and  (29)  indicate  that  the  flow  over  the 
tip  portion  of  the  airfoil  is  of  the  type  referred  to  as  “conical 
flow.”  For  this  type  of  flow  the  values  of  induced  downwash, 
aerodynamic  loading,  etc.,  are  functions  merely  of  the  angle  17. 


831184—49 - 2 
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Busemann  (reference  7),  Stewart  (reference  8),  and  Lager- 
strom  (reference  9)  have  developed  analyses  for  certain  plan 
forms  which  are  postulated  on  the  existence  of  this  type  of 
solution.  In  all  cases  for  which  the  flow  field  is  conical  the 
problem  is  effectively  two-dimensional.  Such  a  simplifica¬ 
tion  reduces  the  analysis  in  this  report  to  a  consideration  of 
a  single  integral  equation  while  in  the  references  just  men¬ 
tioned  complex  variable  theory  can  be  applied  directly. 

Tip  of  swept-forward  lifting  surface. — Consider  the  tip  of 
a  swept-forward  supersonic  lifting  surface  with  uniform 
loading  (fig.  5),  the  angle  8o  between  the  leading  edge  and  the 
X  axis  and  the  angle  5i  between  the  trailing  edge  and  the  x 
axis  both  being  less  than  the  free-stream  Mach  angle  /x. 


Figure  5.— Tip  of  swept  forward  lifting  surface  with  traces  of  Mach  cones,  coordinate  system, 
and  regions  defined  for  equaffons  (32)  and  (34). 

In  carrying  out  the  integrations  it  is  necessary  to  distinguish 
between  the  type  in  which  the  tip  boundaries  are  behind 
the  Mach  cones  and  the  type  in  which  the  tip  boundaries  are 
ahead.  The  analyses  of  these  two  cases  are  of  equivalent 
complexity,  however,  and  can  be  handled  with  equal  facility 
by  the  methods  outlined.  For  all  surfaces  whose  leading 
edges  form  an  apex,  only  the  case  where  the  wing  boundaries 
are  behind  the  Mach  cone  will  be  considered.  A  Cartesian 
coordinate  system  is  chosen  as  shown  so  that  the  origin 


lies  at  the  apex,  the  positive  x  axis  extending  downstream, 
the  y  axis  extending  laterally,  and  the  z  axis  being  directed 
normal  to  the  plane  of  the  plan  form  and  to  the  free-stream 
direction.  The  equations  of  the  sides  of  the  lifting  surface 
are 

y=0 

y=  — X  tan  So=  — ^  x 

and 

y=  —  {x~c)  tan  c) 

The  calculation  of  (p{x,  y,  z)  again  must  be  divided  into 
oases  depending  upon  the  location  of  the  point  P:  (x,  y,  z). 
In  the  residts  listed  below  are  included  the  explicit  expressions 
for  (p  {x,  y,  z) ;  the  induced  velocities,  however,  are  given 
only  in  the  plane  2=0,  as  the  integration  to  obtain  a  general 
e.xpression  is  difficult.  The  velocities  in  the  2=0  plane, 
which  are  sufficient  for  the  purpose  of  this  investigation, 
can  be  obtained  from  a  simpler  integration  since,  for  tbe 
integral  involved, 

lim  r 

2-^0  I{x,y,z)dx= 

This  simplification  was  used  in  the  analysis  of  most  of  the 
lifting  surfaces  investigated.  As  before,  it  is  assumed  that 
the  discontinuity  in  p  is  equal  to  Co.  Moreover,  the  ex¬ 
pressions  for  Wa„o  are  given  in  terms  of  the  variables  17 
and  0)  where 

7;=^^  and  01=1 (30) 
2/ 

In  this  manner  the  solution  is  shown  to  be  conical  in  the 
region  ahead  of  the  I  railing-tip  Mach  cone  (fig.  5).  For 
points  behind  this  Mach  cone  the  flow  is  not  conical  but  a 
function  of  both  r;  and  co. 

Region  7i:  Inside  the  leading-tip  Mach  cone  and  ahead  of 
the  trailing-tip  Alacli  cone.  Integration  of  equation  (17) 
yields  the  result 


J  ^(2',2/,C 


_ ^ 


and,  after  further  calculation. 


xyPx!^  °—pz% 

'  z-y,ix^-^\i/+z^)_ 
(31) 


»“-27rF„J_7  <^vu-i<v<0 

Cod  r  Vl-’t'  1  r  1  j 

27rU,;L  V  h|  + 


(1  +  007?) 

I  (0o+'>7)  L 


,  -1<7,<1  (32) 
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Region  lo'.  Inside  both  tip  Mach  cones.  The  solution  in 
this  region  is  simplified  through  use  of  the  fact  that  for  linear 
differential  equations  any  algebraic  sum  of  solutions  will  be 
another  solution  of  the  equation.  Since  the  differential 
equation  for  the  acceleration  potential  is  linear,  this  property 
can  he  applied  to  obtain  a  solution  for  the  region  I2  by  sub¬ 
tracting  from  the  expressions  given  for  region  /i  correspond¬ 
ing  expressions  in  which  the  variable  5i  replaces  5o  and  {x-c) 
replaces  x.  Thus 


G, 

^  =  2^ 


xy^-x^ 

-arc  tan  ■ — ======= -1- 

z  V?— 

.\2^1 


—arc  tan 


xy 


arc  tan- 


Z^X?  —  (7/^  2^^) 

{x—c)y 


2V(a:— c)^— 2^) 

and  for  — 

[ 


arc  tan 


{x—c)y+{x—cf^—Pz-ei 

(33) 


w 


=  r VW--1.  arc  cosh  n-  + 

^0  It?  I 


■2xFo 


V 


■\^l—da  1  (l-b^oT?)  —  ,  bo 

-  „ —  arc  cosh  ,7^  , — i  — - -  -l-  „  arc  cosh  |-|  — 

00  |(,P0  +  ’?)|  T?  Ol 


arc  cosh 
Qi 


I  (^O  +  T?) 

CO-1-01’7 


Sico-fij  1 


(34) 


lifting  surface  with  trapezoidal  plan  form. — The  linear 
property  of  the  differential  equation  may  be  used  to  advan¬ 
tage  in  determining  the  flow  about  a  trapezoidal  lifting  sur¬ 
face  with  uniform  lift  distribution,  since  the  boundary 
conditions  within  the  plan  form  of  the  airfoil  are  obviously 
satisfied  when  the  acceleration  potential  for  a  triangular 
tip  is  subti'acted  from  the  potential  for  the  rectangular 
surface. 

Suppose  (fig.  6)  the  angle  of  rake  of  the  trapezoid  is  5o  RirI 
that  So  is  less  than  the  Mach  angle  jx.  The  acceleration 
potential  will  be  identical,  over  the  central  portion  of  the 
surface,  to  that  for  the  lifting  surface  of  infinite  aspect  ratio. 
Over  the  parts  of  the  surface  which  are  blanketed  by  the  tip 
Mach  cones  the  flow  will,  however,  be  modified.  Because 
of  symmetry  the  determination  of  this  modifieation  need 
only  be  carried  out  on  one  side  of  the  figure. 

If  the  coordinate  axes  are  chosen  as  shown  in  figure  6,  the 
lateral  boundary  of  the  lifting  surface  is 

,  .  daX 

y=—x  tan  6^=  —  -- 

It  has  been  shown  that  both  the  rectangular  plan  form  and  the 
triangular  plan  form  experience  conical-type  flow  over  the 
region  within  the  tip  Mach  cones.  Thus,  the  variable  t] 
defined  in  equation  (30)  may  be  used. 

Region  Ii :  Inside  the  Mach  cone  originating  at  the  leading- 
edge  tip,  outside  the  Mach  cone  from  the  trailing-edge  tip, 
forward  of  the  trailing-edge  wedge,  and  to  the  left  of  the 
y=0  plane. 


For  —  t<7j<0: 


Wz=0~ 


27rt^n 


Vl  — drti  \ 

Vi  ’71  +  00/ 


^  Co/3  r_E 
27rFo  L  2 


-j-arc  tan  ■■-;+£  arc  cosh  £  — 
VI— >7  00  [j?! 


00  I  (00  +  77)  I J 


(35) 


Uo 


Fujure  <».— Trapezoiclul  lifting  surfaco  with  traces  of  Mach  cones,  coordinate  system,  and 
regions  dofined  for  equation  (35). 


Swept-back  lifting  surface. — As  another  example  of  the 
way  in  which  the  linearity  of  the  differential  equation  may  be 
utilized  to  obtain  further  solutions,  the  induced  vertical 
velocities  for  a  swept-back  wing  will  be  determined  for  the 
(‘ase  in  which  tlie  leading  and  trailing  edges  lie  behind  their 
respective  Mach  cones  (fig.  7).  The  boundaries  of  the  plan 
form  are  given  by  tlie  equations 

00  00 

V=-^jix-c),  y^^jix-c) 

The  flow  will  be  conical  ahead  of  the  trailing-edge  Mach 
cone  where  the  induced  velocities  can  be  expressed  in  terms 
of  the  variable  77.  Behind  the  trailing-edge  Mach  cone  the 
flow  will  not  be  conical  but  will  be  expressible  in  terms  of  the 
c 

variables  and  0=  1  — . 

X 

Consider  first  the  region  of  conical  flow.  In  order  to 
determine  for  a  given  value  of  77  it  is  possible  to  con¬ 
sider  separately  the  induced  effects  produced  by  each  half  of 
the  surface.  But  in  the  region  ahead  of  the  trailing-edge 
Mach  cone,  the  induced  velocities  arising  from  one  half  of 
the  surface  are  given  by  the  formula  for  a  similar  region  on 
the  swept-forward  surface.  For  reasons  of  symmetry  the 
results  for  the  entire  swept-back  lifting  surface  need  only  be 
given  for  values  of  i?  within  the  limits  —  l<t7<^0. 
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Afo 


Figure  7— Swept-back  lifting  surface  with  traces  of  Mach  cones,  coordinate  system,  and 
regions  defined  for  equations  (36)  and  (37). 


In  the  region  where  the  flow  is  not  conical  the  solution  will 
be  built  up  of  a  combination  of  solutions  obtained  from  the 
regions  of  conical  flow. 

Region  h :  Inside  the  leading-edge  Mach  cone,  outside  the 
trailing-edge  Mach  cone,  and  to  the  left  of  the  y=0  plane. 
For  — 1<5;<0 


one  of  technique  in  integration.  The  solutions  for  the  uni¬ 
formly  loaded  sirrfaces,  however,  are  particularly  useful.  By 
methods  of  superposition  these  solutions  can  be  used  to  ob¬ 
tain  the  surface  loading  for  specified  plan  forms  (the  inverse 
problem)  as  will  be  illustrated  in  the  following  section. 


LOAD  DISTRIBUTIONS  ON  FLAT-PLATE  LIFTING  SURFACES  IN  SUPERSONIC 

FLOW 


Infinite  span  wing. — Since  the  vertical  induced  velocity  is 
constant  for  the  supersonic  airfoil  of  infinite  aspect  ratio 
(equation  (22))  and  uniform  load,  it  follows  that  the  airfoil 
is  a  flat  plate.  This  property  distinguishes  the  infinite  aspect 
ratio  problem  from  all  other  plan  forms  considered,  for  the 
load  distribution  must  be  modified  in  the  latter  cases  so  that 
twist  and  camber  are  removed  from  the  wing  to  obtain  a  flat 
plate. 

Denoting  the  angle  of  attack  of  the  airfoil  by  a, 


Moreover, 
and,  setting 


it  follows  that 


Pi—Pu~Pai<Pu—<Pi)  =PoOo 
Pi—Pu_^P 


Ap  2Co 


(38) 


(39) 


Eliminating  Co  between  equations  (38)  and  (39), 


Ap _  4a 


(40) 


ISCgdo  T’  Vl 


J-V 


2  2 
Vi 


m 


dm 


^  Cq^ 
^itVq  _ 


^  arc  cosh  arc  cosh 

PO  I’?  I  ^0 


(36) 

(2-go2-„2)-| 

J 


Region  I2:  Inside  both  Mach  cones  and  to  the  left  of  the 
y=  0  plane.  The  solution  in  this  region  can  be  produced  by 
subtracting  from  the  value  of  Wz=o  given  for  region  C  the 
value  of  Wj.o  given  for  the  same  region  except  that  in  the 
latter  case  5o  is  replaced  by  5i  and  x  by  (x—c).  Thus, 


:  Do/3 


—  arc  cosh 
m 


arc  cosh 


Vi-gi^ 


Although  the  uniformly  loaded  lifting  surface  was  the  only 
prescribed  loading  analyzed,  it  should  be  noted  that  the 
basic  integration  leading  to  a  solution  of  this  type  of  prob¬ 
lem  (equation  (17))  is  in  no  way  restricted  to  a  uniform  load. 
Arbitrary  loadings  that  may  or  may  not  be  analytic  functions 
of  X  and  y  can  be  specified  and  the  problem  therefore  becomes 


The  result  given  in  equation  (40)  is  the  well-known  Ackeret 
expression  developed  in  reference  10.  The  derivation  here 
follows  the  approach  of  Prandtl  (reference  2). 

Rectangular  plan  form. — Since  the  vertical  induced  velocity 
for  the  uniformly  loaded  supersonic  airfoil  of  rectangular 
plan  form  is  not  constant  over  the  portion  of  the  wing 
covered  by  the  tip  Mach  cones,  it  is  necessary  to  modify  the 
load  distribution  within  this  region  in  order  to  get  a  flat  plate. 
The  determination  of  the  required  load  distribution  will  be 
shown  to  depend  on  the  solution  of  an  integral  equation  and 
subsequent  problems  dealing  with  other  plan  forms  will, 
from  a  mathematical  standpoint,  be  similar  in  form. 

The  rectangular  plan  form  will  be  thought  of  as  being  built 
of  superimposed  trapezoidal  lifting  surfaces  with  variable 
angles  of  rake  (fig.  8),  each  trapezoidal  surface  having  a  uni¬ 
form  load  distribution  but  with  loading  allowed  to  vary  with 
the  variable  rake  angle  8. 

Since  the  flow  over  the  part  of  the  airfoil  within  the  Mach 
cone  is  conical,  it  is  possible  to  express  w^^o  as  a  function  of 
t]  where 


Setting 

j3  tan  8  =  8 
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and  using  equation  (35) 

r’''(7'(0)  f-TT-  f’  -^)dd  (41) 

/3  Je=o  \  J-i  Vi  Vi+('/ 

where  C'(6)=^u—<pi  for  the  single  trapezoidal  surface  with 
rake  angle  5. 

M-g 


X 

Figure  8.— Rectangular  plan  form  built  of  superimposed  trapezoidal  lifting  surfaces  with 

variable  rake. 


The  solution  of  the  problem  depends  on  the  determination 
of  a  function  C'  (d)  which,  when  substituted  in  equation  (41), 
will  yield  a  constant  value  of  ^2=0(17) ;  that  is,  a  value  of 
independent  of  the  variable  rj.  Imposing  the  condition  that 

dr, 


the  problem  is  resolved  into  one  of  solving  the  equation 

o=|-  c'ie)de  f” 

dr]Je=o  J-i  Vi  m  +  S 

By  means  of  the  notation 

the  integral  equation  is  written  in  the  form 

o=lim  1";^  C'(e)G,(v,e)ds+^  f'"'  C {e)G,{n,d)de\ 

jj+e  _1 


where  the  singularity  in  the  integrand  necessitates  the  use 
of  the  infinitesimal  «.  The  evaluation  of  the  derivative 
thus  leads  to  the  expression 


'-‘is  [J;"'  f 

C"  ( — 17  —  e)  (’7 ,  —  >7  —  e)  +  C'' ( — ’Z  +  e)  <3'i  (>? )  — +  «)  J 


ViouKE  9.— Hegion  of  integration  showing  line  of  singularity  for  equation  (44). 


It  can  be  shown  from  equation  (35)  that,  if  C'{S)  is  a  con¬ 
tinuous  function, 

lim  f— C"(— >7  — 6)<^i(’7)— ^7— «)  +  C''(~’7  +  «)G’i(’7)— ’7+e)l=0 
e-40  L  -I 


Hence 


0= 


VF=V  r  C'(e)de 


-X 


171  Jo  171  +  ® 
and  the  solution  of  this  equation  is 

a 


C'{d)- 


'  Vs(i— ^) 


(42) 


(43) 


where  61  is  a  constant  to  be  determined  later.  Substituting 
from  equation  (43)  into  equation  (41) 


27rFoW(’7)z=o  1 

f’  rrde 

P  de  r 

/3C.  J 

loV®(l— ®)  J 

VizilL. 

171  (’71  +  ®) 


drii 


(44) 


The  region  of  integration  in  the  7?i®  plane  for  the  double 
integral  of  equation  (44)  is  shown  as  the  cross-hatched  area 
of  figure  9,  a  singularity  in  the  integrand  occurring  along  the 
line  0=— r?i.  Rewriting  the  equation  and  reversing  the 
order  of  integration  in  the  double  integral, 


22L\^Mi-=-2xarcsinV®T-li«^r  dr, 

PGi  Jo  €-4oJ-i  Vi 


r  f'"'  - 

[_Jo  (iji 


de 


+d)My=T) 


r  - 

c/— 77+€  (rji 


de 


+0)V®(1-®) 


] 


The  bracketed  expression  in  this  equation  can  be  shown  to 
vanish  for  all  values  of  r,i  between  zero  and  —1  so  that, 
finally, 


2irVo 


2rr  arc  sin 


—  (7ii8ir 
■  2Va 


(45) 


Since  the  trapezoidal  lifting  surfaces  are  superimposed, 
the  loading  C(d)  over  the  resultant  rectangular  plan  form 
satisfies  the  relation 


dC(e) 

de 


0'(e)  = 


Or 

V®(1-®) 
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Imposing  the  condition  that  (7(0)  =  0  at  0=0,  it  follows  that 

(7(0)  =  2(71  arc  sin  ^[d  (46) 

This  equation  gives  the  incremental  change  of  acceleration 
potential  between  the  upper  and  lower  lifting  surface  of  the 
rectangular  wing.  As  a  result  the  increment  in  pressure  is 

Pi—Pu~Po{‘Pu~<Pi)=2poGi  arc  sin  -yj'd 

Expressing  the  pressure  difference  in  non  dimensional  terms, 

pj^_^4a  Y0  (47) 

IpoVo^  2  n 

The  constant  Ci  may  be  eliminated  between  equations  (45) 
and  (47)  and  as  a  consequence 

A®  8  .  lo 

-^=  — -A-  TT-  arc  sin  ^8 
<Z  Vo  Pv 

Since  the  angle  of  attack  a  of  the  airfoil  is  by  definition 

equal  to  the  final  expression  for  the  loading,  in  coeffi- 

r  0 

cient  form,  over  the  outer  portions  of  the  rectangular  wing  is 


airfoil  must  be  adjusted  properly  to  give  constant  induced 
vertical  velocity.  Superposition  of  trapezoidal  lifting  sur¬ 
faces  with  loadings  varying  with  rake  angle  5  can  again  be 
used  and  the  conical  nature  of  the  flow  employed.  Setting 

n^Pylx 

0=/3  tan  5 

6o=P  tan  00 


equation  (35)  leads  to  the  expression 


2tFowm.^=o=  r*^'c"(0) 
P  Je=Oo 


f  r 

■’>  Vi— dvi  \ 

-1  ’71-F0/ 

(52) 


where  C {d)=(pu—‘Pi  for  the  single  trapezoidal  surface  with 
rake  angle  6. 

The  analysis  in  this  case  follows-  along  lines  directly  analo¬ 
gous  to  that  used  for  the  rectangular  surface.  For  the 
present  configuration  the  loading  function  for  the  superim¬ 
posed  trapezoids  is  given  by  the  relation 


V  (.0- 


(7. 


0o)  (1-0) 


(53) 


Ap 

Y 


8a 


arc  sin 


V0 


(48) 


and  the  integration  to  obtain  can  be  simplified  to 

give,  as  a  final  result. 


The  general  approach  used  to  obtain  this  result  is  similar  to 
that  used  by  Schlichting  (reference  11).  The  error  in 
Schlichting’s  final  result  has  been  noted  by  Busemann  and 
others. 

Lift  coefheient  Cl  for  an  arbitrary  wing  is  defined  by  the 
relation 


ri  —Y—Jl 


Ap 

Q 


dS 


(49) 


PC]  irdd _ _ 

2^ Vo  J e.  (1-0)  ~  2 Fo 


(54) 


The  loading  C{e)  over  the  resultant  trapezoidal  plan  form 
can  be  found  from  the  relation 


dC{d)^__^Ci  _ 

de  ^/(0_eo)'(i-0) 


where 

i  =  total  lift  of  the  wing 
dfS=  element  of  area  on  the  wing 
fS'o=  total  area  of  wing 

For  the  rectangular  wing  the  values  of  Ap/q  over  the  tip 
and  center  sections  are  given  by  equations  (48)  and  (40). 
As  a  result  of  this  integration 


where  A  is  the  aspect  ratio  and  by  definition  equal  to  the 
ratio  of  the  square  of  the  span  and  the  wing  area.  As  a 
final  conclusion  the  lift-curve  slope  of  the  wing  is 


From  the  boundary  condition  that  (7(0)  =  0 
follows  that 


and 


(7(0)  =2(7i  arc  sin 


Pi—Pj,_Ap_iCi 


je—t 


-g-  PoFo^ 


Vi 


0-00 


at  0=00  it 


(55) 


Elimination  of  C\  between  equations  (54)  and  (55)  and 
introduction  of  angle  of  attack  a  for  —  Wj^o/Fo  gives  as 
aerodynamic  loading  over  the  portion  of  the  airfoil  within 
the  tip  Mach  cones  the  expression 


dpL _ 4  / _ 1 _ A 

da~/3V  2PA/ 


(51) 


Ap 

Y 


Set 

— arc  sin 
irP 


/0-0O 

Vl-^o’ 


O<0o<l 


(56) 


Trapezoidal  plan  form.^ — The  results  given  in  equations 
(48)  and  (51)  are  capable  of  generalization  to  the  case  of  the 
flat  plate  having  trapezoidal  plan  form  and  with  rake  angle 
So  less  than  the  Tvlach  angle  of  the  stream.  For  such  a 
configuration  the  airfoil  is  again  blanketed  in  part  by  the 
tip  Alacii  cones  and  the  loading  in  this  outer  section  of  the 


Figure  10  indicates  the  variation  of  the  loading  over  the  tip 
scctioir  of  the  trapezoid.  The  variable  {Pla)(Aplq)  is 
plotted  against  P  tan  5  for  P  tan  So  equal  to  0,  0.3,  and  0.6. 
The  curve  for  P  tan  5o=0  corresponds  to  the  case  of  the 
rectangular  wing  and  shows  results  in  agreement  with  equa¬ 
tion  (48). 
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[B  fan  6 

Fkjukk  lO.~Load  distribution  over  tip  for  various  trapcdzoidnJ  plan  forms. 


By  moans  of  oquation  (49)  together  with  equations  (56) 
and  (40)  the  lift  coofRciont  of  the  trapezoidal  wing  is  expres¬ 
sible  in  the  form 

^  4a/l~2V^'’ 


Introducing  the  aspect  ratio  A  of  the  wing  where 


A= 


c  tan  5o^ 


one  gets  for  lift  coefficient  the  relation 
2a 


Cr. 


From  equation  (57), 


(58) 


da  ^ 


c  c 

1  —  hi  tan  5o— tan  n 


2b 


1  —  t  an  5o 


(59) 


dC  1 

In  figure  11,  /3  is  plotted  as  a  function  of  Afi  for  0o=O, 


and  1.  The  curve  for  6o=0  agrees  with  results  given  by 
equation  (51)  for  the  rectangular  wing.  All  curves  are 


Mo 

T 


FiiH’re  II.  Variiitidii  of  rrduced  lifl-c-urvu  slope  (3  with  reduced  aspect  ratio  pA  for 
various  frapozoidtil  plan  forms. 


terminated  at  values  of  Ai3  for  which  the  lip  Mach  cones 
intersect  on  the  trailing  edge  of  the  wing. 

Triangular  plan  form,  type  1. — The. pressure  distribution 
over  triangular  lifting  surfaces  with  constant  induced  vertical 
velocities  will  be  developed  in  the  following  throe  sections. 
These  plan  forms  arc  indicated  in  figures  12(a),  12(b),  and 
12(c)  and  shall  bo  denoted,  respectively,  as  types  1,  2,  and  3. 
Types  1  and  2  are  actually  special  cases  of  type  3;  namely, 
the  cases  where  one  leading  edge  is  parallel  to  the  free  stream, 
and  where  both  loading  edges  make  equal  angles  with  the 
stream  direction.  Type  3  includes  any  plan  form  which  has 
leading  edges  swept  behind  the  Mach  cone  but  on  opposite 
sides  of  an  axis  drawn  through  the  vertex  of  the  triangle 
and  parallel  to  the  free  stream;  and,  further,  has  a  trailing 
edge  such  that  the  Mach  cones  from  either  tip  do  not  cross 
the  surface  of  the  wing.  The  principal  reason  for  considering 
the  three  types  separately  is  to  show  the  manner  in  which 
the  spanwise  loading  appears  in  the  solution  of  the  problem. 
In  types  1  and  2  the  proper  load  distribution  is  found  readily 
while  the  final  type  requires  a  more  careful  treatment. 

In  order  to  determine  the  load  distribution  over  the  airfoil 
it  will  be  convenient  to  use  a  differential  element  over  which 
the  loading  is  uniform.  The  elements  may  then  be  summed 
and  the  distribution  of  loading  adjusted  so  that  the  induced 
vertical  velocity  at  any  point  on  the  total  lifting  surface  is 
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Ma 


X 

(b) 


Mo 


(0) 


where w{d,ri)  and  ^(0+ AS,?;)  are  the  velocities  induced  by  the 
triangular-tip  surfaces  with  uniform  loading  and  with  tip 
angles  equal  to  8  and  8-\-A8,  respectively.  Applying  a 
limiting  process, 


Aw 
Ae->0  Ay 


=  lim 

A9^0 


wid-\-Ad,ri)~w{9,ri)  bw 
Ad  “'M 


(60) 


It  follows  that  for  the  resultant  lifting  surface  will  be 
evaluated  by  an  integration  with  respect  to  6.  If 
can  be  expressed  in  the  form  of  an  integral  with  respect  to 
r),  the  relation  for  w^^o  will  then  be  similar  to  those  given 
in  equations  (41)  and  (52)  for  the  previous  plan  forms  and 
the  expectation  will  be  that  the  function  C(9)  can  be  deter¬ 
mined  to  give  constant  induced  vertical  velocity. 

The  method  of  attack  just  outlined  is  postulated  on  the 
existence  of  an  integral  expression  for  Such  an 

expi’ession  is,  however,  obtainable  directly  from  the  integrals 
in  equation  (32).  Integrating  these  relations  by  parts 
after  first  differentiating  by  6,  leads  one  to  the  formulas: 

For  — l<r7<0 


bw^^o 

1 

C”  dy)i  "1 

be 

27rFo 

_viv  +  So)  , 

+  Vl— 

and  for  0<5?<1 

bw^^o 

dCo 

rvw  1 

rv  dr)i  ~\ 

be 

2x^0 

[_?7(77-|-i9o) 

Ji  Vl— ’/iM 

Figure  12.— Triangular  flat  plate  lifting  surfaces,  (o)  Type  1,  (b)  Type  2.  (c)  Type  3. 

constant.  For  the  triangular  plan  forms  it  is  possible  to 
assume  that  conical  flow  exists  and  the  analysis  may  be 
carried  out  using  the  angular  coordinates  that  have  already 
been  introduced. 

Figure  13  shows  the  elemental  lifting  surface  to  be  used. 
The  sides  of  the  element  extend  back  from  the  tip  of  the 
Mach  cone;  making  angles  8  and  5-f-A5  with  the  positive 
X  axis  or  free-stroam  direction.  Corresponding  to  previous 
notation,  the  relations  6=^  tan  8  and  0-}-A0=/3  tan  (6-1-A6) 
are  used.  The  vertical  velocity  induced  by  the  element  of 
surface  may  be  denoted  by  Aw  and  it  follows  that 


If  the  elements  are  summed  over  the  type  1  triangular 
wing,  induced  vertical  velocity  is 


Vi  —  _i_  T’  _ _ 

viv  +  Q)  J-l  ’7i^(’?i  +  0)  Vl— ijil 


dd 


(63) 


where  C{d)=(pu—ipi  for  the  element  at  5=arc  tan  If 
Wz=o  is  constant,  then 

by)  ^ 


Aw—w{d-{-Ad,  rj)  —w{d,  ij) 


and  from  this  criterion  the  function  C{d)  will  be  determined. 
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Thus,  using  methods  similar  to  those  introduced  in  the 
development  of  equation  (42), 


1  C{e)de 
V  Jo  V+6  ’JVo  V  +  6 


^  C(e)de 
Jo  V  +  6 


(64) 


The  general  solution  of  this  equation  is 


C{6)  = 


G,e+C2 

V0(0o-0) 


where  Ci  and  C2  are  constants.  Since,  however,  the  Kutta- 
Joukowski  condition  requires  that  loading  vanish  along  the 
edge  0=0,  it  follows  that  (72  =  0  and  the  required  loading 
takes  the  form 

=  (65) 


If  equation  (65)  is  substituted  into  equation  (63),  vertical 
induced  velocity  can  be  calculated  from  the  expression 

2TrVoWz^a _ — If  r®® _ ddd _ I 

"  -n  jo  (r)  +  0)V0(eo-0)"^ 

P"  ede  __  p  d77i _ 

Jo  V®(®0— 0)  J-1 ’?l^(’?l  +  0)  Vl  — 

The  region  of  integration  in  the  iji,  0  plane  for  the  double 
integral  of  equation  (66)  is  shown  in  figure  14  for  the  case  in 
which  —  A  singularity  in  the  integrand  of  the 

double  integral  exists  along  the  line  0-l-i7i  =  O.  Reversing 
the  order  of  integration,  equation  (66)  may  be  rewritten  as 

^  Jo  {d-{-Tl)  ^l6{6o - 6)  J-l  7?lVl  - Vi^ 


p»  d0 

p  drf 

p«  d0 

Jo  (0-b77i)  V0(0o-0) 

J-e(i77i-v/l  — 771^,, 

lo  (0-h7,.)V0(0O-0) 

(67) 


The  single  integral  in  equation  (67)  has  a  singularity  at 
0— — 7?  since  — 0o<7'>?^O  and  ri  therefore  lies  inside  the  region 


0 


Figure  14.— Region  of  integration  showing  line  of  singularity  for  equations  (66)  and  (78). 


of  integration.  A  corresponding  singularity  occurs  in  the 
second  of  the  double  integrals  at  0=  — 771.  Consider,  there¬ 
fore,  the  integral 


j: 


d0 _ r  p 

(0+v)  ^0(00-0)  €->0  LJo 


d0 


(0-1- 77)  V^O — 0 


p»  do  1  .gg, 

J— ’?+«  (0-\-7j)  V0(0O 0)  — 


The  indefinite  integral  is 


^  — 7700  T^^0~l~  2770 — 2V  (  —  7700 — 77^)0(00—^ 

■\I—r)0a—Tf  0-1-77 


so  that  the  definite  integral  is 


lim  ~i= 
.->0  y- 


-7700- 


In 


gp  I 

-0, 


In  [— 27700  — 600  — 277^— 2677— 2-v/(— 7700— 77^)(  — 77  — f)(0o-l-77-|-6)]  ^  ) 

[—2)70o-|-e0o— 277^-1-2677— 2  V(— 7700— 77^)(—77^6)(0o-l-77— 6]  (—6)  ) 


The  value  of  this  expression  is  0  and  equation  (67)  therefore 
becomes 


27rVo'M77-o^  _  f”®”  ..  ^^>71  r””..  _  (09\ 

J-l  771  \/i  — 77i^)o  (04-771)  —  0) 

Since,  in  this  region  of  integration,  —  l<y77<  — 0o  it  follows 
that 


j; 


d0 


(0+7?,)V0(0O-0) 
1 


—  77 100  4- 000-1- 27710 

'7;,0«4'77.^  ^''"2V(47p+7,A)0(0„-0) 


"Do 

Jo 


\  77  (00-1- 77) 


and 


drf 

°  2F0J-I  771-v'l— 77iV7)i(7;i4-0o) 


(70) 


The  integral  of  equation  (70)  can  be  transformed  by  means 
of  purely  algebraic  substitutions  into  a  form  that  integrates 
immediately  into  complete  elliptic  integrals  of  the  first  and 
second  kind.  However,  in  the  consideration  of  the  type  3 
plan  form  it  will  be  necessary  to  resort  to  other  methods  of 
transformation,  so  that  a  more  uniform  approach,  employing 
Jacobian  elliptic  functions,  will  be  used  throughout.  (See 
reference  12.) 

The  quartic  under  the  radical  in  equation  (70)  is  first  re¬ 
duced  to  an  expression  of  the  type  appearing  in  elliptic  in¬ 
tegrals  of  canonical  form.  This  is  accomplished  by  succes¬ 
sive  application  of  the  transformations 


771 


k+ls  ,  k 
—  and  s=-r 

14-s  t 


where  k  and  I  are  chosen  so  as  to  destroy  the  odd  powers  of 
the  variable.  By  means  of  these  transformations,  induced 
velocity  becomes 
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pq  (i-F)F 

2Vo  ^k{do—k){l-k^)ji  (i_Pi2)Y(l-«2)(Pi2_i) 


where 


e„ 


The  integration  of  equation  (71)  will  be  performed  after 
first  considering  two  parts  such  that  Wz^o='Wi-{-W2  where 


^ _  {i—V)k‘^  r k _ _ 

'^^~2Vo  ■\/~k{eo—k){i—P)ji 


where  E'  is  the  complete  elliptic  integral  of  the  second  kind 
with  modulus  k'=-\ll — k^.  Equation  (72)  can  be  further 
simplified  by  writing  k  in  terms  of  so  that 


_  /3(7i  /2(l+Vl-eo^)  r., 

2Fo  V  00^ 


where  the  modulus  of  E'  is  and  k= 


l-Vl-go^ 


For  the  loading  in  question 

Pi—Pu  =  Pa{‘Pu  —  ipi)  —  poGi 


,  /  e 

'  V 


_  0Ci  (1— A:^)P 


W2=  — 


2Fo  ^k{e^-k){\-B)}i  {l-kV)^{l~f){kH-^-\) 


This  separation  is  prompted  by  the  fact  that  the  integral  for 
Wi  is  expressible  in  terms  of  elementary  functions  after  the 
simple  transformation  s.  The  results  of  such  an  integra¬ 
tion  lead  to  a  value  that  is  zero  at  the  lower  limit  and  infinite 
at  the  upper  limit.  However,  an  inspection  of  the  original 
integral  in  equation  (70)  shows  that  is  finite  so  the 
infinity  obtained  for  Wi  must  be  canceled  by  a  corresponding 
infinity  of  equal  magnitude  in  w^.  The  actual  proof  of  this 
statement  necessitates,  of  course,  treating  the  combined 
expressions  as  an  indeterminate  form  where  the  upper  limits 

of  the  integrals  for  Wi  and  are  replaced  by  y-t-  e  and  the 

k 

limit  is  taken  as  e  approaches  zero. 

Introduce  now  in  the  integration  of  W2  Jacobian  elliptic 
functions  and  set 


so  that 


t=sn{u,  k)  =snu 


dt=cnu  dnu  du 


The  expression  for  becomes 
{l—B)kH 

^Vo  ^k{do—k){l~k^) 


J’K+iK' 
/< 


where  K  and  K'  are  the  complete  elliptic  integrals  of  the 
first  kind  with  respective  moduli  k  and  k'=-sjl — k'^.  Inte¬ 
grating  and  combining  with  W\,  one  has 


{\-k‘^)kH  f  1 
2Vo  ^k{e,-k)(l-h)  W 


1  r  .  j^snuenu  .1)*+'^*'' 


where  E(u)  is  the  incomplete  elliptic  integral  of  the  second 
kind.  After  substitution  of  the  limits,  induced  vertical 
velocity  is 


Wz=o=  — 


2Fo  ^k(e,-k)(i-k^} 


Pl—Pu_^P 
^PoVo^  2 


_2q  /_J 

~Vo^  V^o- 


By  means  of  equation  (73)  the  constant  q  may  be  eliminated 
and 

2a  ,  /20(1-VTF^ 

q  ^E'M  d,-e  ^  ^ 


O  ■£  .4  .6  .8  LO 

13  fan  6 

Figure  15.— Load  distribution  over  triangular  plan  forms  of  type  1. 
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Figure  15  shows  the  variation  of  ^  ~  with  /3  tan  8  for  values 

of  jS  tan  5o  equal  to  0.3,  0.6,  and  0.9. 

From  equations  (74)  and  (49)  the  lift  coefficient  of  the  right 
triangle  wing  with  trailing  edge  normal  to  the  free-stream 
direction  can  be  determined.  It  follows  that 

—  Vl  —  ^0^) 

Since  the  aspect  ratio  A  is  given  by  the  equation  A=26ol^, 
equation  (75)  can  be  used  to  find  the  lift-curve  slope  as  a 

dC 

function  of  A.  In  figure  16  a  plot  of  /3  as  a  function  of 
^A  is  given. 


2^Vo  f-e,  L^(j+^  + J-i  + 

or,  since  C{6)  =  C{~-d), 

w(,),=o=  - \j¥=^  +  J_i 

(76) 

The  function  0(6)  in  equation  (76)  must  give  a  constant 
value  for  w(ri)i=o  so  that  will  vanish.  Imposing 

this  condition  it  can  be  shown  that  a  solution  is  given  by 
the  relation 

C, 


and,  after  substitution  in  equation  (76) , 


(77) 


Mo 


ttVoWz^o r°°  yi — Tj^do  _ 

“dC'i  “  Jo 

C""  _  do  _  p  dr)i 

Jo  ojdo^  —  sO-i  Vl  — ’?i^) 


(78) 


The  integration  of  equation  (78)  is  to  bo  performed  under 
the  assumption  that  — (9o<C»?<CO  so  that  the  region  of 
integration  in  the  r/i,0  plane  for  the  double  integral  is  as 
shown  in  figure  14.  Reversing  the  order  of  integration  in 
the  double  integral,  equation  (78)  may  be  written  in  the 
form 


0A 


Figure  Hi.— Variation  of  reduced  lift-curve  slope  0  with  reduced  aspect  ratio  BA  for 

da 

plan  forms  of  typo  1. 


Triangular  plan  form,  type  2.— Figure  12(b)  shows  the 
symmetrical  type  of  triangular  plan  form  considered  in  this 
section.  The  semivertex  angle  is  8o  and  is  defined  by  the 
relation 


TrVoW.^o 

■ 


VT 


'X‘ 


dd 


0  (0^-7,^)  V^o^-V^ 
r-®"  p _ dd _ 

J-i  Vl— ’ll"  Jo 


C”  P" _ 

J-eoT/iVl — Jo  {8^  — 


dd 

-ni^)  V<9o^— 


Evaluation  of  integrals  of  the  form 

dd 


1= 


('eo 

■Jo 


is  accomplished  by  means  of  the  substitution 

8 

do 

After  substitution,  the  integral  becomes 


0o=/3  tan  8o 

The  loading  element  used  in  the  previous  section  can  be  used 
again  and  equations  (61)  and  (62)  are  applicable  directly. 
Because  of  the  symmetry  of  the  figure,  it  is  necessary  merely 
to  insure  the  constancy  of  w^^o  over  the  left  half  of  the  wing 
in  order  that  the  entire  wing  be  a  flat  plate. 

Summing  the  elements  over  the  type  2  triangular  wing, 
induced  vertical  velocity  over  the  portion  of  the  wing  for 
which  —  1<C’7<C0  is 


and,  by  straightforward  integration, 


J= 


0  for^,<l 


(79) 
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This  result  shows  that  the  second  double  integral  of  equation 
(79)  vanishes  (since  Oa'^vi)  as  does  also  the  single  integral 
in  the  equation.  For  the  remaining  double  integral,  however, 
6o<C.'n\  and 


Setting 


fiCi  f  ^ _ 

1 

2= - 

Vi 


(80) 


equation  (71)  transforms  to 


pCi  _ 


Introducing  the  modulus  k  =  do  and  making  the  substitution 
z=sn(u,  k)=snu 
one  gets  the  expression 


W, 


''2Vjc^ 

pc, 


K+iR’ 

K 


2VJc- 


E' 


where  the  prime  again  refers  to  the  complementary  modulus 
k'  =  of  the  complete  elliptic  integral.  Since  k  =  d(, 

pc. 


where  A;'  =  Vl  — 

For  the  loading  in  question 

Pi—Pu_^P_ 


2Foeo' 


E' 


(81) 


2C, 


PoVo^ 


FoS/0o^-e^ 


so  that,  eliminating  C,  between  this  equation  and  equation 
(81), 

Ap  _  iadp^ 


a  ^^Oo^—e^E' 
P  Ap  . 


(82) 


Figure  17  shows  the  variation  oi  —  —  with  p  tan  S  for  values 

of  p  tan  So  equal  to  0.3,  0.6,  and  0.9. 

From  equations  (82)  and  (49)  it  is  possible  to  find  the 
expression  for  lift  coefficient  of  a  triangular  or  delta  wing. 
Thus: 

2ird„a 


(83) 


Since  aspect  ratio  of  the  wing  is 


Mo 


where  the  modulus  of  E'  is  k 


This  result 

6 


agrees  with  that  obtained  in  another  manner  by  Stewart 

dC 

(reference  8).  In  figure  18  a  j^lot  is  given  of  /3  ^  ^ 


function  of  pA. 

Triangular  plan  form,  type  3. — Figure  12  (c)  shows  the 
plan  form  now  to  be  considered.  Relative  to  the  x  axis  or 
free-stream  direction  the  sides  of  the  triangle  form  the  angles 
6o  and  6i  so  that  the  total  vertex  angle  is  5o+Si  =  2A.  The 
variables  0o  and  6,  are  also  introduced  satisfying  the  relations 
0o=P  tan  5o,  d,  =  P  tan  5,.  The  same  loading  element  that 
was  used  for  type  1  and  type  2  triangles  may  be  used  and 
equations  (61)  and  (62)  apply.  It  will  then  be  necessary  to 
determine  the  distribution  of  load  so  that  the  induced  vertical 
velocity  over  the  plan  form  is  a  constant.  Since  this  induced 
velocity  must  be  the  same  on  both  sides  of  the  5=0  axis, 
two  equations  result: 


lift  coefficient  becomes 


A= 


4^0 

P 


C^-- 


TaA 


(84) 


For 


'ro(’jh-o= 


1<’7<0 

p  r  rvi 


2-a-Vc 


r  r. 

J—6,  J-i 


drii 


-1  7;i^(7?i  +  0)  Vl  ■ 


C(e)d0 


-Vi  -I 


(85) 
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The  evaluation  of  and  Hi{di,do)  is  accomplished 

in  the  same  manner  as  has  been  used  previously:  first,  a 
reduction  of  the  quartic  under  the  radical  to  canonical  form 
and  second,  transformation  by  means  of  Jacobian  elliptic 
functions  lollowed  by  direct  integration.  Since  the  cacu- 
lations  for  both  equations  are  quite  similar,  only  in  the  case 
of  do)  will  the  details  be  mentioned. 


0  .a  I.e  2.4  3.2  4.0 

HA 

Figure  is.— Variation  of  reduced  lift-curve  slope  fi  '-j''  with  reduced  aspect  ratio  fiA  for 

(la 

plan  forms  of  typo  2. 


and  for  0<ij<l 


__dv\ 


-ni 


^  C{e)d9 


(86) 


From  the  solutions  to  the  problems  of  type  1  and  type  2  it  is 
possible  to  construct  a  solution  of  the  more  general  problem 
by  expressing  the  loading  function  in  the  form 


Ad-\-B 


(87) 


where  A  and  B  are  constants  that  can  be  determined  in  terms 
of  1^2=0  from  equations  (85)  and  (86).  Equation  (87),  in 
conjunction  with  equations  (85)  and  (86),  yields  the  expres¬ 
sions 

'W2_o=2^  [AHi{do,  di)ABH2{do,9i)\ 


where 


Wi^o — 2^  [  do)  A BH2{6i,  do) 


drti 

Vi-y/ A  — vi'‘){vi  —  0i){vi  +  So) 

dyi 

V  (1  —  ’/i^)  (»!i  —  ^i)  (’ll  +  ^o) 


(S8) 

(89) 

(90) 


i-dod,-^l{i-do^){i-9,^) 

®  00  +  ^1 

^_i-eoer+V(i-eo^)(i-ei^) 

^0+^1 

The  introduction  of  the  symbols  R  and  k  defined  as 

b~a  _ 

~  V (1—0^)  {do— a)  (01 +  a) 

,  i-dob 

^^'To-b 


(‘quation  (89)  reduces  to 


1 


— h)  —  l]dt 

(T-a¥) 


(91) 

(92) 

(93) 

(94) 

(95) 


The  integrand  divides  naturally  into  two  parts,  one  containing 
even  powers  and  one  containing  odd  powers  of  t  in  the  nu¬ 
merator.  The  laftei'part  integrates  into  elementary  functions 
after  substituting  t=u^  and  equation  (95)  thereby  becomes 


F-Jldt _ -| 

{l-a^f)^/{t^-l){l-kH^)J  ^  ^ 


Setting 


■j:  (.- 


(f—\)dt 


aH^)-,/{f-l){l-kH^) 


and  substituting 


x=sn{u,  k)  =snu 


one  gets 


1-aWuJ 


du 


If  sn  7=^'  h  DOW  may  be  written  as 

^  .  dny  rK+iK'  k^snyenydnysrdu  , 

^  k^snyeny  Jk  l  —  k^sTdysrdu  ^ 

or  _ 

h=-K'-i^^^2mu,y)W^'^'  (97) 

where  11  (tt,^)  is  the  fundamental  elliptic  integral  of  the 
third  kind. 
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The  evaluation  of  11  {u,  k)  is  best  achieved  by  means  of  its 
expression  in  terms  of  theta  functions  and  zeta  functions. 
Thus 


nKT)=2  log 


Q{u—y) 

e(li+7) 


]ruZ{i) 


and  the  bracketed  term  in  equation  (88)  is 


Q{K-\-iK' —y)Q(K-\-y)  , 
e  {K-  7)  e  {K+  iK'+'yy 


iK'Ziy) 


The  theta  functions  are  quasi-periodic,  that  is,  they  satisfy 
the  relations 

Q(u+2K)==e(u) 


From  this  property,  together  with  the  fact  that  ©(«)  is  an 
even  function,  it  follows  that 

e{K+iK'-  y)e{K+y)  _ 

GiK-y)QiK+iK'+y) 


Moreover,  since 

Z(y)=Eiy)—y  ^ 


there  results 


UiK+iK',y) -U(K, 7)  =iK'E{y) -iy  i  ^ 

ir+.y|Fi4|i+K'[£(7)-7|]|  (98) 

The  expression  for  equation  (96)  can  now  be  written 

Hi=a^Rk^ | ['E'(7)-7 | -a~RkK 

(99) 

where  the  moduli  are  k  for  the  nonprimed  functions  and 
^'=-y/l  — F  for  the  primed  functions.  By  defijiition,  7=arc 

sn  fQ,  where  F  is  the  incomplete  elliptic  integral  of 

the  first  kind  with  argument  ^  and  modulus  k. 

In  the  same  notation,  the  equation  for  is  as  follows: 

a\l-e)K- (1  -«  J  it'  j  r ^ 

.M) 


and 


H,= 


kRb 


E{y)—y 


K 


Ki-0]| 


Formulas  (99)  and  (100)  can  now  be  combined  with  equa¬ 
tions  (79)  to  give 


A- 


and 


Vo'^^z^O  (a  _a\  I 

I3E'~  1^0  Veo+01 


(101) 


2G^ 

+  ^i 


fhere 


G= 


1  —  V  (1  —  ^0^)  (1  ^1^) 


(102) 

(103) 


and  E'  is  the  complete  elliptic  integral  of  the  second  kind 

with  modulus  yX  —  Q^. 

From  equations  (87),  (101),  and  (102) 


Ti—Vu 


ApJ2C{e)^2 

2  V 0^  V(^i"l“^)  (^0 — 


2a 

pE 


/  2G  ^ 

(6o — 6i)  04-200^1 

\  60+ 6, 

_V(^i  +  ^)  (^0  6)_ 

(104) 


It  should  bo  remarked  that  the  slope  of  the  loading  curve 
is  zero  at  6=0.  Figures  19  (a),  19  (I)),  and  19  (c)  show  the 

variation  of  —  ~  with  B  tan  5  for  values  of  /3  tan  6i  equal  to 
a  g 

0.3,  0.6,  and  0.9,  respectively,  and  for  /3  tan  8o  equal  to  0, 
0.3,  0.6,  and  0.9. 
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From  equations  (49)  and  (104)  the  lift  coefficient  for  a 
type  3  plan  form  is  obtainable.  Two  cases  will  be  developed 
here:  first,  when  the  trailing  edge  of  the  wing  is  perpendicular 
to  the  stream  direction;  second,  when  the  trailing  edge  of 
the  wing  is  perpendicular  to  the  line  of  symmetry.  The 
first  configuration  may  be  referred  to  as  a  skewed  wing  while 
the  second  configuration  may  be  referred  to  as  a  symmetrical 
delta  wing  at  an  angle  of  sideslip.  Thus  for  a  skewed  wing 

Ci=-g7^  -y/ {do-\-6i)20  (105) 

where  G  is  given  by  equation  (103)  and  E'  has  the  modulus 
-\ll  —  G^.  This  result  agrees  with  that  given  by  K.  C. 
Roberts  in  an  abstract  in  reference  13. 

For  the  more  practical  case  of  the  delta  wing  at  an  angle 
of  sideslip,  figure  12  (c),  the  lift  coefficient  can  be  expressed  as 

^  2onr  .  jG  tan  A 
<?!,=  ^-r  cos  Ay - - - 


where  A  is  the  angle  of  sideslip  and  2A  the  angle  between  the 
leading  edges,  and  0  is  expressed  in  terms  of  0o  and  di  which 
are,  in  turn,  expressed  in  terms  of  A  and  A  by  the  following 
equations 

6o=/3  tan  (A-|-A)  ) 

6/i=i3  tan  (A-A)  I 

Since  the  pressure  distribution  has  been  computed  only  for 
wings  with  leading  edges  behind  the  Mach  cone  springing 
from  the  apex  and  with  a  trailing  edge  ahead  of  the  Mach 
cones  from  the  wing  tips,  formula  (106)  is  valid  only  for 
cases  where 

m+A<90°n 

AH-A<m  V  (108) 

A— A>0  ) 


(106) 


These  restrictions  are  practically  always  met,  however,  for 
angles  of  sideslip  likely  to  be  encountered  in  flight. 
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shown  as  a  function  of  sideslip  and  A.  The  figure  shows 

dC 

that  up  to  15°  of  sideslip  remains  practically  constant. 


90° 

a  fx 
4  -A^  O 


A,  degrees 

Figure  20.— Variation  of  reduced  lift-curve  slope  /3  ^  with  angle  of  sideslip  A  for  some 
plan  form.s  of  type  .3. 


Equation  (106)  is  plotted  in  figure  20  where 


dC^ 

da 


is 
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National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  April  I4,  1947. 

REFERENCES 

1.  Prandtl,  L.:  General  Considerations  on  the  Flow  of  Compressible 

Fluids.  NACA  TM  No.  805,  1936. 

2.  Prandtl,  L.:  Theorie  der  fiugzeugtragflugel  in  zusammendriick- 

baren  medium.  Luftfahrtforschung,  Bd.  13,  Oct.  1936. 

3.  Bateman,  H.:  Partial  Differential  Equations  of  Mathematical 

Physics.  Dover  Publications  (New  York),  1944. 

4.  Webster,  Arthur  Gordan:  Partial  Differential  Equations  of  Mathe¬ 

matical  Physics.  G.  E.  Stechert  &  Co.  (New  York),  1933. 

.5.  Hadamard,  J.:  Lectures  on  Cauch3'’s  Problem  in  Linear  Partial 
Differential  Equations.  Yale  University  Press,  1928. 

6.  Lamb,  Horace:  Hydrodynamics.  Dover  Publications  (New 

York),  1945. 

7.  Buseman,  A.:  Infinitesimal  Conical  Supersonic  Flow.  NACA 

TM  No.  1100,  1947. 

8.  Stewart,  H.  J.  The  Lift  of  a  Delta  Wing  at  Supersonic  Speeds. 

Quarterly  of  Applied  Mathematics,  Vol.  IV,  No.  3,  Oct.  1946, 
pp.  246-254. 

9.  Lagerstrom,  P.  A.:  The  Application  of  Analytic  Extension  in  the 

Solution  of  Problems  in  Supersonic  Conical  Flows.  Abst.  No. 
7,  pub.  No.  3,  JPL,  GALCIT,  1946. 

10.  Ackert,  J.:  Air  Forces  on  Airfoils  Moving  Faster  than  Sound. 

NACA  TM  No.  317,  1925. 

11.  Schlichting,  H.:  Airfoil  Theory  at  Supersonic  Speed.  NACA  TM 

No.  897,  1939. 

12.  Whittaker,  E.  T.,  and  Watson,  G.  N.:  A  Course  of  Modern 

Analysis.  4th  Edition,  Cambridge  University  Press,  1940. 

13.  Roberts,  R.  C.:  On  the  Lift  of  a  Triangular  Wing  at  Supersonic 

Speeds.  Abst.  No.  382,  Bulletin  of  the  American  Mathematical 
Society,  Vol.  52,  No.  11,  pt.  1,  Nov.  1946. 


U.  S.  GOVERNMENT  PRINTING  OFFICE:  1S49 


Positive  directions  of  axes  and  angles  (forces  and  moments)  are  shown  by  arrows 


Axis 

Force 
(parallel 
to  axis) 
symbol 

Moment  about  axis 

Angle 

Velocities 

Designation 

Sym¬ 

bol 

Designation 

Sym¬ 

bol 

Positive 

direction 

Designa¬ 

tion 

Sym¬ 

bol 

Linear 
(compo¬ 
nent  along 
axis) 

Angular 

Longitudinal . 

X 

X 

Rolling 

L 

Y—-*Z 

Roll 

u 

D 

Tifttfirfll 

Y 

Y 

M 

Z - ¥X 

Normal  . 

Z 

Z 

N 

X - *Y 

Absolute  coeflficients  of  moment 
n _  ^  ^ 

(roUmg)  (pitchmg) 


O 

(yawing) 


Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5.  (Indicate  surface  by  proper  subscript.) 


4.  PROPELLER  SYMBOLS 


D  Diameter 

p  Geometric  pitch 

pID  Pitch  ratio 

V'  Inflow  velocity 

V,  Slipstream  velocity 

T 

T  Thrust,  absolute  coefficient  Gt— — stw 

pn^D* 

Q  Torque,  absolute  coefficient 


P 


O, 


■n 

n 

$ 


Power,  absolute  coefficient  Cp= 


P 

pn^D^ 


Speed-power  coefficient = 


^jpV^ 

yPn^ 


Efficiency 

Revolutions  per  second,  rps 


Effective  helix  angle =tan~^ 


( 


2-Krn/ 


1  hp=76.04  kg-m/s=550  ft-lb/sec 
1  metric  horsepower =0.9863  hp 
1  mph=0.4470  mps 
1  mps=2.2369  mph 


5.  NUMERICAL  RELATIONS 


1  lb=0.4536  kg 
1  kg=2.2046  lb 
1  ini=  1,609.35  m=5,280  ft 
1  m=3.2808  ft 
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